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We determine the regime where the widespread classical field description for quantum Bose gases 
is quantitatively accurate in Id, 2d, and 3d by a careful study of the ideal gas limit. Numerical 
benchmarking in ID shows that the ideal gas results carry over unchanged into the weakly interacting 
gas. The optimum high energy cutoff is in general shown to depend strongly on the observable in 
question (e.g. energy, density fluctuations, phase coherence length, condensate fraction). This 
explains the wide spread of past results. A consistent classical field representation with less than 
10% deviation in all typical observables can be given for systems at temperatures below 0.0064 
degeneracy temperature in Id, and 0.49 critical temperature in 3d. Surprisingly, this is not possible 
for the 2d ideal gas even at zero temperature because mean density, density fluctuations and energy 
cannot be simultaneously matched to the quantum results. 


I. INTRODUCTION 


The quantum mechanics of a wide variety of physical 
systems can be quite accurately described by an appro¬ 
priately chosen ensemble of complex fields (also called, 
classical or c-fields) Idl . [^ . Examples include quan¬ 

tum gases of ultracold atoms, coherent light fields, and 
solid state polariton systems. A common feature is the 
appearance of collective behavior such as high amplitude 
phase fluctuations and superfluid defects, that strongly 
fluctuate away from the mean held. Though the term 
“classical” is used, we are talking about the opposite 
regime to the usual gas of classical particles. Here it is 
the collective field that has classical properties such that 
each member of the ensemble could be non-destructively 
tracked, while the particles lose their individual identity. 
Examples of such a ppr oaches include classical held en¬ 
sembles l4l|, the stochastic Gross-Pitaevski 

Muation [3, Izl, [8[, the truncated Wigner representation 
[ 9 I-HIII for ultracold atoms, and the open stochastic clas¬ 
sical held equations for polaritons [l2, El] • Related ap¬ 
proaches for fermions include stochastic mean held the¬ 
ory in e.g. heavy ion collisions E3, El] , and effective held 
theories for the pairing order parameter . 

In the absence of sufficient in situ experimental reso¬ 
lution, the approach is also commonly used like a hight 
recorder to give information on the dynamics of the sys¬ 
tem before its detection in destructive time-of-hight im¬ 
ages. Its applications are growing in importance given 
advances in the experimental investigation of sponta¬ 
neous superhuid defects and phase huctuations, such as 
d, [H-lll, Isl - fil . Idl . 1^ . IHl . In quantum many-body sys¬ 
tems with collective nonlinear phenomena, such ensem¬ 
bles of complex helds are often the only practical way to 
obtain theoretical information on huctuations, full distri¬ 
bution functions, and - especially - on typical single real- 



However, their use has usually been accompanied by 
lingering doubt on whether the results are quantitative or 
qualitative. From an operational perspective, two major 
contributing factors to that have been (1) a visible de¬ 


pendence of some results on the high energy cutoff that 
is chosen, and (2) different optimum cutoff values found 


the literature El, E3, El E3, Ei, EH, m, H il 
l69j| . The aim of this paper is to identify a regime where 
c-helds are in fact quantitatively accurate, so that they 
can be used there with conhdence. 

Qualitatively, the condition for the applicability of 
classical helds to bosons is that the relevant physics 
can be captured by considering only the highly occupied 
sing le-particle modes, without the need for a condensate 
[6^ . Poorly occupied modes are not described well, and 
those above an energy cutoh need to be discarded to 
avoid pathological behavior such as the UV catastrophe 
known since late 19th century physics. 

However, the matter of just where to draw the line 
and how accurate the description is, has been a matter of 
much contention and ambiguity. The history of applying 
classical helds to ultra cold atomic gases teaches us that 
accuracy has depended quite strongly on the choice of 
the high energy cutoff and the observables studied. Past 
numerical benchmarking E ^ E^H ^P- careful 

comparisons to experiment TZ, Um|, and also analyti¬ 
cal pH . [ 2 ^ and purely mathematical studies [IHl of var¬ 
ious single observables have found that it is possible to 
achieve good to very good agreement but the details of 
the recipe vary from study to study [l|, E3, E3, EH, EH]. 

Here, we intend to clarify these dependencies, and will 
show that under the right conditions the classical held 
approximation can be treated as more than just a quali¬ 
tative guide, but gives predictions that are correct within 
small error bounds for a wide range of observables. 


We will concentrate hrst on the case of an ideal gas 
as a baseline, reasoning that well described interacting 
regimes can be found at temperatures that are already 
well-described in the ideal gas. Then we will conhrm that 
accuracy seen in the ideal gas carries over into the weakly 
interacting regime under appropriate conditions. We will 
work in the local density approximation (LDA) in the 
thermodynamic limit. That is, we will consider pieces 
of the gas cloud having a certain local density, which 
allows us to remain general in terms of trap geometry. 
In the LDA, it is natural to work in the grand canonical 
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ensemble (GCE), where the rest of the system acts as the 
particle and thermal reservoir. Such a model underpins 
more general behavior, and it will be seen that several 
important conclusions can be reached. 

In Secini we will describe our approach. Further, in 
Sec lIII Al we will find the temperature dependent “eigen” 
cutoffs that allow the classical fields to correctly match 
the density and one other observable. Subsequently, in 
Sections IIII Bl and IIVI we will determine the resulting er¬ 
rors in other observables and the cutoff that minimizes 
the systematic error across the whole range of observ¬ 
ables. This will tell us about the temperature range over 
which an accurate complex field description of the sys¬ 
tem is possible. Finally, via numerical benchmarking of a 
weakly interacting gas in \d to the exact Yang and Yang 
solution [ 3 ^ we will show in Sec IVIl that the ideal gas re¬ 
sults carry over largely unchanged into that regime. We 
conclude in Sec IVIIl 

II. APPROACH 
A. Classical field description 

The essence of the classical fields method is to re¬ 
place annihilation (creation) operators au (aj)) of single 
particle modes in the field operator by complex ampli¬ 
tudes which is warranted when occupation is 

macroscopic. Then we can write: 

’J'W = ^a/cV'fe(x) I I (1) 

k I fcgc J 

where ipk{^) is the wave function for the fcth mode and 
C denotes the low energy subspace. Since we will be 
considering uniform sections of the gas, plane wave modes 
fc = k are the most convenient, with momentum cutoff 
kc so that only modes |k| < kc are included in C. 

In general, it should be understood that ^(x) corre¬ 
sponds to an ensemble {...} of complex field realizations, 
each with its own set of amplitudes ■ The full ensemble 
preserves the gauge symmetry of the quantum thermal 
state that corresponds to a set of many experimental re¬ 
alizations. This is despite the ’’virtual” symmetry break¬ 
ing done by each member of the ensemble similarly to a 
single experimental realization . 

B. Parameters 

The properties of the uniform dilute gas can be en¬ 
capsulated by two dimensionless parameters. The first is 
7 = with density n and contact interaction strength 
g, and the second is a reduced temperature r which 
depends on the density n, but not on the interaction 
strength. We choose the thermal de Broglie wavelength 

At = J as our length scale, so that the reduced 


temperature is 

Td 2 ti h? ^ ’ 

Here, Td is the usual quantum degeneracy temperature in 
d dimensions that corresponds to one particle per region 
of volume A^. It is a natural scale for our investigation 
because then r = I corresponds to the point at which the 
highest mode occupation is 0 ( 1 ), and this constitutes the 
intuitive ultimate upper bound on temperature for which 
classical field descriptions make sense. 

It is convenient to also scale the cutoff in these units: 

/c = fcc (3) 

A value of /c = 1 corresponds to a cutoff at the plane 
waves with thermal de Broglie wavelength At. We will 
henceforth work in the following units: At = 1 and H = 
m = 1, where m is the mass of particles. Note that the 
cutoff in terms of single particle energy is 

ec = 7rkBT f^. (4) 

In the ideal gas limit (7 0) that we consider first, 

there is only one physical parameter characterizing the 
system - the density-dependent reduced temperature r, 
and one technical parameter fc for the classical fields 
description. Phase space density equal to one occurs 
at T = tt) = {1.539,1.443,1.368} in Id, 2d and 3d re¬ 
spectively, while the BEG critical temperature in 3d is 
T = Tc = 0.5272. 


C. Observables 

The great majority of experiments concentrate on low 
order observables such as phase, density or their fluctu¬ 
ations. We will analyze the following: 

1. n - density. 

2 . e - kinetic energy per particle. 

3. Ipg ” phase grain length. 

This is the size of a coherent region, which we will cal¬ 
culate via Ipg := i / dz ('k'l'(O)^(z)) = Jdzg^^^(z). 
In the quasicondensate regime, when g^^^z) ~ 

Ipg equals the phase coherence length 

4. g*'^^(0) - normalized local density fluctuations. 
While these are of much theoretical interest, they are 
rarely measured in situ because imaging resolution is 
usually much worse than the intrinsic density correlation 
length of the system. 

5. uc - coarse-grained density fluctuations. 

This quantity is defined as uq '■= vaTN/{N) = 
n J dz (z) — l] -1-1, where N is the atom number in a 
region much larger than the density correlation length. In 
contrast to g^^^(O), this intensive thermodynamic quan¬ 
tity often appears in experimental work [ 2 ^ and it 
gives the ratio of the measured fluctuations in a pixel 
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to Poisson shot noise. It is equal to the static structure 
factor at |k| = 0, i.e. S'(O) . 

6. po - condensate fraction. 

7. Qr - coherence half width. 

In the presence of a true condensate, Ipg (and uq) ceases 
to be a good thermodynamic quantity, diverging because 
—>■ oo) = Po- In light of this we need another 
measure of the width of phase fluctuations, and will 
define it by the half width of the peak of i.e. 

g(i)(a^) = 1(1 + 

It is worth noting that the kinetic energy per particle 
in itself is not a typical subject of measurement, but its 
consideration has here its own justihcations. If typical 
observable quantities are described correctly, but e is not, 
then this will quickly come out as errors in the dynamics. 


D. Ensemble 

A major consideration in our work here has been to 
remain independent of trap geometry. This basically re¬ 
quires working in the local density approximation (LDA). 
As an example of variations with geometry that can occur 
without an LDA approach, optimal energy cutoffs found 
on the basis of the distribution of condensate fraction for 
a whole cloud in the canonical ensemble were 0.29fcBT in 
a uniform box, but l.OfesT for a harmonically trapped 
gas [2l|. In the end in Sec. IVK. we will see that the re¬ 
sults of the LDA approach can be largely reconciled with 
the harmonically trapped canonical ensemble results. 

When considering a relatively uniform section of a 
larger gas, it is not only possible, but also essential to 
work in the grand canonical ensemble (GCE) rather than 
the canonical one. In such a situation the rest of the 
system acts as a particle and thermal reservoir, while 
the uniform GCE section describes the properties that 
are local to the region. This approximation is acceptable 
provided the physical length scales such as Ipg are shorter 
than the length scale on which the density changes. Such 
conditions generally prevail for quasicondensates or a 3d 
gas above the condensation temperature. 

Use of the GCE in a truly condensed system such as 
the 3d gas below or the finite-size 2d gas at extreme 
low temperatures, requires some care and background to 
get our bearings. It is known that for the ideal gas the 
usual thermodynamic equivalence between ensembles is 
lost in the presence of condensation. Particularly glaring 
differences are seen in the fluctuations of condensate frac¬ 
tion between the canonical and grand canonical ensemble 
— a matter that has been much studied [ziMlzilTl 
and is sometimes known as the “fluctuation catastrophe” 
for the GCE. In fact, a uniform condensed system in the 
GCE has anomalous fluctuations of the number of con¬ 
densed particles (i.e. their variance grows faster than the 
mean number), which implies that some quantities such 
as UQ diverge. Technically this signals the point of the 
breakdown of the theory [79l - [^ , but in reality this kind 
of behavior cannot actually occur. Physically the growth 


of diverging quantities is braked by other effects. Usu¬ 
ally, the causes can be traced to either a breakdown of 
the thermodynamic limit due to hnite size effects, or a 
suppression of fluctuations due to interactions (see [1^ 
for a detailed discussion). 

The primary difference between the grand canonical 
and microcanonical or canonical treatments of an ideal 
condensed system has been pointed out quite early Q 
by studying the ground state number fluctuations. They 
are huge in the GCE (the occupation Nq of the ground 
state is exponentially distributed P{Nq) ~ 
small in the other thermodynamical ensembles. In con¬ 
trast, there is no such difference for excited level occu¬ 
pations. This suggests that the majority of observables 
are not pathological. Even the mean condensate frac¬ 
tion does not diverge nor break equivalence between en¬ 
sembles, unlike its fluctuations. Hence, it is legitimate 
to benchmark classical fields in the GCE provided that 
we exclude from consideration those observables that are 
known to be deviant. In particular, when condensation 
is present neither condensate fraction po nor the main 
coherence decay described by ar are pathological, so we 
will use these instead of uq and Ipg. 

The suppression of anomalous condensate fluctuations 
due to interactions can occur even at very weak interac¬ 
tions. This can be seen from a simple argument: Con¬ 
sider the GCE partition function of the condensate mode: 

OO 

Zo{p,T)= e^t^No-CogNl)/T (- 5 ^ 

Nq—Q 

with g the interaction strength, and Cq a geometry- 
dependent factor [^. In this form it is now a Gaus¬ 
sian distribution of the condensate occupation with mean 
Nq = p/2gCo and relative condensate number fluctu¬ 
ations ANo/No = y/T/{2CogN§) = ^T/{2E,^t)- It 
means that the relative magnitude of the number fluc¬ 
tuations is related to the ratio of the temperature to 
the interaction energy of the entire system, The 

latter very quickly suppresses the grand canonical fluc¬ 
tuation catastrophe as the size of the system becomes 
appreciable, leaving only a tiny low temperature region 
at r < 1/Eint with anomalous fluctuations, that shrinks 
as T —0. 

The above considerations are distinct from the sepa¬ 
rate matter of what ensemble should be considered for 
the entire system. If one were to nondestructively fol¬ 
low a single realization of the system over time and as¬ 
sume ergodic evolution, then the correct ensemble would 
be the microcanonical one that has the system isolated 
from particle and energy exc hang e. This has been con¬ 
sidered in many works [ 3 , [H, iTol. ITSL [7§| . On the other 
hand, actual experimental studies usually deal with an 
ensemble over many independent realizations created by 
cooling a new cloud each time, and independently mea¬ 
suring each destructively. Then the fluctuations of the 
number of particles between different realizations can in 
fact be of the same order as the mean number of parti- 
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cles over the whole experimental series. Due to the large 
number fluctuations between shots, a sequence of single 
clouds is likely to be more closely described by the GCE 
than the CE. 

To wrap up this section, once correlation length scales 
are short enough for the LDA to be valid, the approach 
used here is relevant also to global properties of the sys¬ 
tem when the GCE fluctuation catastrophe is suppressed. 
This can happen because of any of the following: (1) lack 
of a true condensate, (2) observables that do not depend 
on fluctuations of the condensate fraction, (3) weak but 
sufficient interaction for Eint of the entire system to be 
large compared to the temperature, (4) an experimental 
data set that consists of an ensemble of many indepen¬ 
dent realizations, except for the cases with strong post¬ 
selection on particle number. Such conditions prevail in 
a very wide range of systems of interest. Keeping these 
in mind, let us proceed. 


E. Benchmarking 

We will compare the classical field predictions for the 
observables in Sec. Ill C| to the exact Bose gas values in the 
thermodynamic limit. For the ideal gas they can mostly 
be obtained analytically. 

To proceed, the LDA approximation requires hrst the 
density to be correct, in the sense that an ideal gas with 
density n (i.e. reduced temperature r) should be com¬ 
pared to a classical field ensemble with the same den¬ 
sity. This is also essential in practice regardless of the 
LDA, since n is the most basic observed quantity in ex¬ 
periments. To match ideal gas and classical field den¬ 
sities, firstly chemical potentials and respec¬ 
tively, must be chosen. A sum over Gibbs factors gives 
the exact Bose gas density and the density es¬ 
timate , fc) in classical fields as functions of 

their grand canonical chemical potentials. We invert 
these, and with the help of Eq. m obtain 
and Other observables e.g. and 

can then be expressed as functions of r and 

fc as well. 

In general, for the Bose ideal gas, r and the choice 
of units specify all properties of the system. In classical 
fields, in addition to r, the system description requires 
a technical parameter fc. There, we can fit both densi¬ 
ties and to r, but also we can make one other 
quantity agree exactly by an appropriate choice of fc- 


III. OBSERVABLE-DEPENDENT ACCURACY 

A. Single observable “eigen” cutoffs 

Fig. m shows how such cutoffs matched to different ob¬ 
servables (which we will call eigen cutoffs) behave as a 
function of temperature. 


The density is already matched due to the LDA as ex¬ 
plained above, and is not shown. We have also not shown 
results for (0) because it is always correctly predicted 
to be g^^\0) = 2 for every cutoff in the ideal gas. This 
property will not hold any more when interactions are 
present. Indeed, then the local density fluctuations man¬ 
ifest a dependence on cutoff. 

The high temperature behavior is qualitatively simi¬ 
lar in all dimensions. The eigen cutoffs matched to en¬ 
ergy per particle ff and to coherence half width /“'■ rise 
to constant values, while the eigen cutoff matched for 
density fluctuation drops to zero (this will be com¬ 
mented on later in Sec. 1^3). The fc^” takes intermediate 
values and is almost constant. An unexpected feature is 
the similar behavior of cutoffs corresponding to and e 
rather than the and Ipg that are more related physi¬ 
cally. 

The crossover to low temperature behavior is around 
T = 1, as expected. In the low temperature 

regime, most eigen cutoffs collapse to a common value 
(0.436 and 0.564, in Id and 2d, respectively), except for 
ff which prefers the higher values 0.653 and 0.724. In 
3d, the cutoffs at r —>■ 0 are 0.783, 0.753, and 0.653 for 
e, Ur, and po, respectively. 

Below critical temperature in 3d, the eigen cutoff for 
condensate fraction has a constant value. This comes 
about because the critical temperature in classical fields 
is cutoff dependent, = [4/c]“^/^, while in the Bose 
gas it is Tc = [C(3/2)]“^/^ = 0.5272 with C(3/2) the Zeta 
function. The condensate fractions are directly related as 
Po‘^^ = [1 - (t/tc-)^] and p^'o^'^ = [1 - ("r/Tc^^)^]. Hence, 
fPo — _ 0.65309 makes po’^'^ and p^f^^ equal for all 




FIG. 1. (Color online) Matched eigen cutoffs fc for several 
observables as a function of temperature t {e - blue line, Ur 
- purple line, Ipg - green line, uq - red line, po - orange line). 
The top panels (a), (b) show Id and 2d cases, respectively, 
and the bottom panels present the 3d situation with (d) a 
magnification of the critical region. The Bose gas critical 
temperature tc is marked as a vertical dashed line, while the 
black solid line shows the fc value below which condensation 
of classical fields occurs. 
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FIG. 2. (Color online) Variation of the relative errors da 
of observables with cutoff /c at representative high and low 
temperatures r. Colors like in Fig. [T] Top row: Id, second 
row: 2d, third row: 3d above rc, last row: 3d below rc- 


T <TC- 

Two other noteworthy points are that: (*) in 2d, the 
eigen cutoff fc = l/v^ that gives the correct phase 
grain length Ipg does not depend on temperature, and 
(**) the wave-like behavior of in Id (as well as in 2 d), 
that comes from oscillations of 5 ^ 1 ) (z) with distance, is 
caused by the sharp cutoff in momentum space in classi¬ 
cal fields. 


B. Relative errors of single observables 


Now, how does a non-optimal choice of fc affect the 
observables, and their systematic error? This is very rel¬ 
evant for practical considerations. For one thing, in a 
nonuniform system, when the cutoff is matched in one 
spatial region, it is good to know the sensitivity of results 
in other regions with a different density on this choice 
of fc- Furthermore, we need this information to judge 
how good the classical fields are in describing the system 
overall. 

The relative error Sa of an observable a is: 


Sa{T,fc) 


Aa 

a 


( a^''^Hr,fc) _ \ 

y Q;('d)(T) J 


( 6 ) 


Its cutoff dependence is shown in Fig. [21 The first obser¬ 
vation is that the relative error of energy per particle has 


an opposite trend to the other quantities. The resulting 
mismatch turns out to be the strongest restriction on the 
fc range for which all Sa errors are small. 

Secondly, in Id the known fact [ 2 ^ that and 

do not depend on cutoffs at low r, is reflected 
in small errors in Ipg, uq, and a,.. However these errors 
are no longer small in higher dimensions. As temper¬ 
ature drops, the Sa{T, fc) except for Sp^, collapse onto 
curves that stay invariant with r and remain steep (the 
T = 0.08 and r = 0.05 panels in Fig. |2|). In other 
words, observables remain sensitive to cutoff all the way 
down to zero temperature in 2d and 3d. 


IV. GLOBAL ACCURACY 

What does it take to match all, or at least to be close to 
all typical observables? Let us consider the global error 
estimator 


RMSa,p,...{T, fc) = Y -I- { 5 ^ + ... 

This is a root mean square of the relative errors of chosen 
observables a, (3, etc. Each relative error will, by defini¬ 
tion, be less than RMS. The main aim of the function 
RAIS will be to catch inaccuracy in any observable. 

We have studied the RAISa,p,... with all the observ¬ 
ables that we have been considering. Moreover, we also 
took various combinations of them. It turns out that 
when we include just uq and e, all relevant features that 
were seen with larger sets of observables are covered. 
This happens because these quantities are the most “ex¬ 
treme” in terms of the behavior of eigen fc and of the 
values and trends of 5a. This is seen in Figs. [DEI Also, 
the pair (e, ug) includes observables of 2nd and 4th or¬ 
der in 41 , which are the two main classes measured in 
experiments. We will use them to define the quantity: 


RMS{t, fc) = (7) 

that will be our indicator of the overall accuracy and 
applicability of the classical fields approximation. Below 
Tc in 3d, the condensate fraction po will be used instead 
of Uq. 

Minimizing Eq. ([7]) at a given temperature will give the 
optimal cutoff momentum and minimum error indicator 
minRMS. For example, a minRMS value below 0.1 (i.e. 
< 10% error in observables) is often satisfactory and we 
will take it as a guideline. 

Fig. El shows the results for the Id gas. Global error 
RMS is very large above the degeneracy temperature 
T = 1. For low temperatures it falls to zero, as one 
would hope. According to our 10% guideline, classical 
fields give acceptable results up to r = 0.0064. The best 
choice of fc is fairly invariant with temperature in this 
region, being in the range (0.649 ± 0.043). In fact, if we 
choose the average value of fc, we will be close to absolute 
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FIG. 3. (Color online) Summary results for the Id gas. The 
left panel shows the dependence of the global error estimator 
RMS, based on e and uq, on cutoff fc for several values of 
T = {0.008 brown, 0.08 orange, td gray, 8 purple}. The 
top right panel shows the minimal value of RMS achieved at 
the optimal cutoff shown in the lower right panel. An addi¬ 
tional dashed branch indicates a less optimal local minimum 
of RMS. 


minRMS regardless of temperature or density. At high 
T an extra second branch appears that is associated with 
a local minimum of RMS with large errors in e and small 
in UG- It is not of practical importance for us because it 
is less optimal. 

Fig.|4]shows the results for the 2d gas. The behavior at 
low temperature is surprisingly unfavorable. RMS never 
falls below 0.333. This is a consequence of an inability 
to satisfy both observables uq and e. Their relative er¬ 
rors (t,/ c) become stuck on the curves shown in 

the fourth plot of Fig. [^whenever r < 0.08 and do not 
cross near zero error. One wonders whether this situa¬ 
tion {minRMS well above 10% as r —>■ 0) is repeated 
for other different sets of observables? It turns out that 
even the pairs {e,lpg) or (e, a^) will lead to similar large 
minRMS values. In fact, no combination that includes 
e and any other observable will work well, because the 
Sa{T, fc) curves are invariant. The crucial and a priori 
not so obvious conclusion is that in 2d, in the small tem¬ 
perature, ideal gas regime the classical fields description 
gives at best only a qualitative description of the gas, 



T 


FIG. 5. (Color online) Summary results for the 3d gas. The 
solid lines correspond to RMSe.uc in the temperature region 
T > Tc and dot-dashed lines to RMS^.p^ in the region r < 
Tc- The left panel shows the dependence of the global error 
estimators on cutoff fc'. t = td gray, r = 0.85 brown, r = 
0.55 purple, r = 0.45 red, r = 0.25 orange, and r = 0.05 
pink. The black vertical line indicates the critical cutoff for 
TC. The right panels are as in Fig. [31 with the Bose gas 
critical temperature tc marked with a vertical dashed line 
and solutions below tc as dot-dashed lines. 


and a description that is quantitatively correct across ob¬ 
servables is unreachable. The matter of whether this is 
alleviated once interactions become important warrants 
further study. 

Fig-Elshows the results for the 3d gas. The area above 
critical temperature behaves analogously to low dimen¬ 
sions. However, around the critical temperature, the 
RMS curve narrows and the accuracy of classical fields 
becomes very sensitive to the choice of the cutoff fc. This 
is related to the fluctuations uq growing to infinity at tc- 
As such, it may be related to the inequivalence of the con¬ 
densed ideal gas ensembles and may be an effect that is 
readily removed by finite size or interaction effects. In 
the condensed regime below tq, the RMS curve widens 
out again while classical fields rapidly become accurate 
with RMS < 10% below r = 0.486 rc. 


V. DISCUSSION 



FIG. 4. (Color online) Summary results for the 2d gas. De¬ 
scription as in Fig. (3) The r = 0.008 and r = 0.08 lines in 
the left panel overlap. 


Several points can be addressed on the basis of the 
ideal gas results, before considering an interacting gas. 


A. Nonuniform gases 

So far, we have been fully focused on the local density 
approach here in order to obtain results that are appli¬ 
cable for general inhomogeneous cloud geometries. 

A very convenient aspect of what we have found is that 
the best cutoff value opt/c is practically constant in the 
whole region where classical fields are a good description 
(say, minRMS <10%). This can be seen in Figures [3] 
andO where in this region, opt/c G (0.645,0.653) in Id 
and opt/c G (0.778,0.783) in 3d. The best low tempera- 






















7 


ture cutoff in 2d is also a constant opt/c = 0.639 — see 
Fig. m Even beyond this best region, the opt/c value is 
almost constant until values of r ~ 1 are reached. For a 
nonuniform gas at a temperature T, the reduced temper¬ 
ature scales with density as t(x) oc l/[n(x)]^/'^. Those 
aspects ensure that the optimum cutoff for all sections 
of the gas is practically the same, provided only that the 
bulk of the gas is quantum degenerate (i.e. r(x) < 1). If 
it isn’t, then the description is not accurate anyway. Op¬ 
erationally, this all means that the best cutoff to choose 
regardless of the density profile of the gas is: 

kc = \/2TnnkBT ; £c = 7r(opt/c)^fcB7’- (8) 

So either we take the low temperature cutoff or it doesn’t 
matter anyway. 

The case of the uniform GCE is in fact quite well 
matched to the trapped canonical ensemble (CE) gas 
that was mentioned in III D1 despite the apparently dif¬ 
ferent framework of the problem. This is because the 
dominant central bulk of such a trapped gas is effec¬ 
tively a uniform open system in the LDA. The cutoffs 
found for the harmonically trapped canonical ensemble 
of the ideal gas based on condensate fraction distribu¬ 
tion at low T correspond in our notation to val¬ 
ues of = {0.56,0.72,0.84} in 1, 2, and 3d, re¬ 

spectively. These are quite close to the GCE values of 
opt/i'"^^^ = {0.65,0.64,0.78} found here (Figs. [51[n]). 
This further reinforces the view that results obtained 
with the LDA are also relevant for nonuniform gases, 
even when the entire cloud does not have particle ex¬ 
change with an environment. 

A certain exception is the canonical ensemble in a box 
whose cutoffs were also studied in na and found to be 
much lower= {0.30,0.47,0.65}. This indicates 
that this is a special case which describes very different 
physics. The matter of which ensemble should be used 
to describe the recently achieved box potentials 
is still open. If the interaction is not too strong, then 
shot-to-shot fluctuations in energy and particle number 
can be appreciable and so a grand canonical approach 
may be warranted for the whole gas (if one is concerned 
with ensemble rather than time-averaged properties). 

B. Breakdown mechanism at high temperatures 

The reason for the drop of to zero at high t pro¬ 
vides an instructive example of how the classical field de¬ 
scription breaks down. Generally the explanation comes 
down to different statistics of particle numbers TVk in the 
modes. Both the fully quantum Bose gas and the classi¬ 
cal field have an exponentially decaying particle number 
distribution: 

P(Afk) oc (9) 

in each mode. However, in the exact treatment, A^k can 
only take on discrete values 0, 1, ..., while in the classi¬ 


cal fields the non-integer part of the distribution is also 
needed. This peculiarity very strongly increases fluctua¬ 
tions especially when the bulk of the distribution is in this 
region, i.e. the mean number of particles is N-^. < 0(1). 

The above observations are transfered to uq in the fol¬ 
lowing way: For the exact Bose gas, the distribution is 
Poissonian when A^k "C 1, giving var[A^^‘'^^] = (A^^'^^) 
for each mode, while the exponential distribution in clas¬ 
sical fields gives var[7V^^'^^^] = Due to having 

independent modes, uq = /X]k(-^k), and 

in the exact treatment ug —>■ 1 directly. To obtain the 

fcf) 

same with classical fields, occupations A}} ~ 1 are nec¬ 
essary to make var[A^k] ~ These are much greater 

than in the Bose gas. So to also simultaneously match 
overall density of the many-mode gas, the cutoff must be 
made much lower than the Bose gas momentum width 
2ttIAt so as to get the same area under the distribution 
of density in k-space. From ([2|), this immediately implies 
<C 1. With such a great modification of A^(k), cor¬ 
rectly matching additional observables like e with classi¬ 
cal fields becomes out of the question. 

A similar breakdown can be expected whenever the 
physics is captured by low-occupied independent modes. 
For example, such discrepancies were seen between ex¬ 
periment and classical fields in the quantum Bogoliubov 
regime of the interacting gas at very low temperatures 


VI. CROSSOVER TO THE INTERACTING GAS 

An obvious question is whether the ideal gas results 
carry over into the interacting gas. To address this, we 
have benchmarked the classical field description in Id sys¬ 
tem with the Yang- Yang exact solution for the uniform 
interacting Bose gas for a sequence of increasing in¬ 
teraction strengths that cover the crossover from the ideal 
gas to an interaction-induced quasicondensate. 


A. Procedure 

The exact values for n = N/L, as well as system en¬ 
ergy E in a segment of length L can be obtained via the 
self-consistent numerical solution of the inte gral equa¬ 
tions given in the original Yang-Yang paper [30|. The 
Hellmann-Feynman theorem was used by Kheruntsyan et 
al. to obtain 5^^^(0) = — (dP/dg)^ j. from the Yang- 
Yang solution for pressure P [H, [s^ , which can be read¬ 
ily evaluated numerically. For the contact-interacting gas 
the expression for the interaction energy in the system 
is Pint = \gn'^Lgk^\Q). From this, one obtains the ki¬ 
netic energy per particle: Ekin = {E — Eint)/N. The 
coarse-grained density fluctuations can also be found via 
UG = idn/dg)rp, based on the expression for var[V] 
in Hi]. 






To obtain classical field results, we generate ensem¬ 
bles of classical field realizations d>(a;) using a Metropo¬ 
lis algorithm, in a way conceptually similar to the work 
of Witkowska et al. [1^ but using grand canonical en¬ 
semble weights g[i^N{’i’)-E{'i’)]/kBT^ numerical lat¬ 

tice is chosen to have a box of length L with periodic 
boundary conditions that is wide enough for the density 
and phase correlations to decay to zero before wrapping 
around. The number of points was 2^*^, which is easily 
sufficient for the maximum numerical lattice wavevector 
to be many times larger than the cutoffs kc imposed on 
the field in k-space. This ensures that no aliasing prob¬ 
lems appear for the evaluation of the interaction energy 
term, as has been discussed in the context of the padded 
lattice in the PGPE and truncated Wigner methods [s^. 
Classical field values for observables at a given cutoff are 
calculated using 10 ^ ensemble members. 

For each cutoff, the observables are benchmarked 
against exact Yang-Yang values for systems having the 
same values of T, 5 , and density n as the classical field 
ensemble. It remains true for the interacting gas that 
UQ and Ekin have the most extreme behavior among the 
set of observables that now also include the interaction 
energy per particle and The latter two have a 

cutoff-dependent behavior that is somewhat similar to 
uc- Hence, we continue to use the same global accuracy 
indicator m as for the ideal gas, using the kinetic energy 
per particle Ekin and coarse-grained fluctuations ug- 
minRMS and opt/c are obtained by fitting a func¬ 
tion to the cutoff-dependent values of RMS{t, fc) at a 
given r and 7 . We use the square root of a parabola 
because it is a good candidate for describing the fc- 
dependent behavior of RMS{t, fc) near the minimum. It 
marries the approximately linear behavior of and Sue 
in this region that is seen in FigEl with the expression 
© for RMS. We use data from an fc range of about 
±0.05 around the minimum. Frror bars are obtained 
by splitting the field samples into Ns smaller subensem¬ 
bles, calculating subensemble values of minRMS^'"^ and 
opt/i*^ in the same way for each, and invoking the cen¬ 
tral limit theorem to estimate the uncertainty in the full- 

ensemble values to be Aopt/^ = \J var[opt/i*^]/As and 
AminRM S = •\/var[mmi?M 5 '(d]/iV 5 . 


B. Results 

We have carried out the above benchmarking for the 
reduced temperature r = 0.00159 and a range of inter¬ 
action strengths 7 = 5 /n from 2 x I0“® to 0.005 in the 
dilute interacting gas. These are experimentally realistic 
parameters. The local bunching g^^^(O) changes over this 
range from 1.976 in the very weakly interacting limit to 
1.02 at 7 = 0.005. This indicates that we move from an 
almost perfect ideal gas on the left deep into the strong 
quasicondensate regime on the right, where almost all ef¬ 
fects are dominated by the interaction mean field. The 



Y 

FIG. 6 . (Color online) A preview of the situation in the Id 
interacting gas. Here, r = 0.00159, and y = gjn increases 
to the right. Top panel: the change in uq, obtained from 
the Yang-Yang lSOll exact solution, compared to the ideal gas 
7 — >■ 0 value (red). Middle panel: minimal value of RMS as 
in Figs. 13151 with la statistical error bars from an ensemble of 
10 “^ samples, and the ideal gas value shown as the horizontal 
line. Bottom panel: corresponding optimum cutoff fc and 
its ideal gas value. One sees that while the observable ua 
changes by two orders of magnitude, the ideal gas values for 
cutoff and accuracy carry over onto the interacting gas. 


coarse-grained density fluctuations change by two orders 
of magnitude over this range, as plotted in the top panel 
of Fig. El 

The results of this foray into the interacting gas are 
shown in the other panels of Fig. El The ideal gas values 
for cutoff carry over onto the interacting gas unchanged, 
to within available statistical precision. The global ac¬ 
curacy minRMS actually improves. One concludes then 
that in this regime at least the optimum cutoff and de¬ 
gree of accuracy found in the ideal gas applies very well 
to a wide swath of the interacting gas as well. This is not 
an a priori obvious result, but certainly a convenient and 
encouraging one for those who want to make calculations 
using classical fields. 


VII. CONCLUSIONS 

To conclude, we have judged the goodness of classi¬ 
cal fields for describing the ideal Bose gas in Id, 2d, 
and 3d using all the usually measured observables. We 
have shown that 10 % or better accuracy for the whole 
set of observables simultaneously is possible in Id up to 
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temperatures of T = 0.0064 with the cutoff prescrip¬ 
tion fee Ri 0.65(;|^) and in 3d up to T = 0.49with 
kc Ri 0.78 {■^)- The essence of the matter can be cap¬ 
tured by the indicator RMS based on kinetic energy per 
particle and coarse grained density fluctuations, which 
are the observables that are the hardest to mutually sat¬ 
isfy. 

In 2d, we have found a surprising feature that classi¬ 
cal fields remain incapable of properly describing all the 
observables together in the ideal gas even as T —>■ 0. 
One suspects that finite size effects and/or weak inter¬ 
actions may improve agreement here. The indication is 
that something is going on in 2d that warrants further 
study. 

When a system is correctly described with a classical 
ensemble of complex fields as here, the observation of 
many “intrinsically quantum” effects that rely on wave- 
particle duality or a discretization of the basis is ruled 
out. This includes things such as stronger-than-classical 
correlations, Heisenberg uncertainty relations, mode en¬ 
tanglement, EPR and Bell inequality violation, anti¬ 
bunching, and noncommuting observables. All in line 
with the difference between classical optics on the one 
hand and quantum optics and quantum information the¬ 
ory on the other. Thus, for parameters in which the 
weakly interacting Bose gas is described by the classical 
field to some level of RMS, observation of the above in¬ 
trinsically quantum effects with typical observables will 
also be suppressed to a level of the same order as RMS. 


Of course, large RMS is not sufheient to imply quantum 
effects. 

Two results lead to optimistic conclusions for the prac¬ 
tical application of classical fields to ultra cold gases. 
Firstly, the optimum cutoffs in the ideal gas are almost 
unchanged with r in the whole region where accuracy 
is good. This means that even for a nonuniform cloud 
with a common global temperature, a single cutoff value 
is close to optimal in the entire degenerate region. This 
goes a long way towards pacifying one of the leading prac¬ 
tical worries. Secondly, our study of the crossover into 
the interacting gas in Sec. IVII shows that the cutoff that 
optimizes the ideal gas is also valid for a part of the inter¬ 
acting gas, including a region where the quasicondensate 
is dominated by interactions. The degree of accuracy 
seen in the ideal gas is also preserved. This is a nontriv¬ 
ial but very encouraging result. A more detailed study of 
the situation for the whole range of interaction strengths 
in Id is in progress and will be reported on in future. 
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